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Abstract 

We calculate the resonant electronic Raman scattering for the Falicov-Kimball model near the Mott transition on 
a hypercubic lattice. The solution is exact, and employs dynamical mean field theory. 
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The Falicov-Kimball model [1] has two kinds of par- 
ticles: conduction electrons, which are mobile and lo- 
calized electrons which are immobile. The Hamiltonian 
is (at half filling) 
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where c] (ci) creates (destroys) a conduction electron 
at site i, fj [fi) creates (destroys) a localized electron 
at site i, U is the on-site Coulomb interaction between 
the electrons, and t* is the hopping integral [2] (which 
we use as our energy unit). The symbol d is the spa- 
tial dimension, and (ij) denotes a sum over all nearest 
neighbor pairs (we work on a hypercubic lattice). The 
model is exactly solvable with dynamical mean field 
theory [3] when d ^ oo (see [4] for a review). 

The formalism for calculating the Raman response 
was originally developed by Shastry and Shraiman [5]. 
The expression for the Raman response is 
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for the scattering of electrons by optical photons (the 
repeated indices a and /3 are summed over). Here 
e''°' denotes the polarization vector of the incident 
and outgoing photon, £j(j^) refer to the eigenstates 
describing the "electronic matter", and Z is the par- 
tition function. The coupling of the photon to the 
electronic system is treated with the linear coupling 
of the vector potential to the current operator j^ = 
^^ 9e(k)/9kc,C|j^Ck, and the quadratic coupling of the 
vector potential to the stress-tensor operator 7^/3 = 
^^^^'^ e{]f.) / d]f.ad]i.j3cl.C]^, with e(k) the bandstructure 
and Ck the destruction operator for an electron with 
momentum k. The scattering operator then becomes 
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with the sum I over intermediate states. 

We have evaluated these expressions for the Stokes 
Raman response, with an incident photon frequency 
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uji , an outgoing photon frequency coo , and a transferee! 
photon frequency Q — u>i — ujo- The procedure is com- 
phcated, and involves first computing the response 
functions on the imaginary time axis, then Fourier 
transforming to imaginary frequencies, and finally per- 
forming an analytic continuation to the real axis [6]. 
Our calculations include effects from nonresonant di- 
agrams, from resonant diagrams, and from so-called 
mixed diagrams. We analyze three different symme- 
tries of the incident and outgoing light: Aig- with 
both incident and outgoing polarizers aligned along 
the diagonal e^ — e° = (1, 1, ...), _B]^„ with the inci- 
dent light polarized along one diagonal e* = (l,l,...), 
and the outgoing light along another diagonal e° = 
(1, —1, 1, —1, ...), and B2p- with the incident light po- 
larized along e* — (1, 0, 1, 0, 1, 0, ...) and the outgoing 
light polarized along e° = (0, 1, 0, 1, 0, 1, ...)• The Ai„ 
response has contributions from all types of processes, 
the B-^g. response is nonresonant or resonant, and the 
-B2g response is purely resonant. 

The Falicov-Kimball model on a hypercubic lattice 
has a Mott transition at half filling when U = \/2. 
The transition is to a pseudogap-like phase, because 
the infinite tails of the noninteracting density of states 
(DOS) (which is a Gaussian) do not allow the system 
to have a true gap of finite width. Instead the DOS is 
equal to zero only at the chemical potential, and there 
is exponentially small DOS within a "gap region". We 
examine the system just on the insulating side of the 
Mott transition at U — 1.5. 
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Fig. 1. The (Stokes) resonant Raman response for three differ- 
ent symmetries as a function of O (at T — 0.05 and U — 1.5; 
just on the insulating side of the Mott transition). The curves 
are for different values of the incident photon frequency. 

Our results for the Raman response as a function of 
the transfered frequency il appear in Fig. 1. Note how 
there are strong resonant effects, including the sharp 
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Fig. 2. Resonant profile for Cl — 0.5 and various cJi > 0.5. 
The curves correspond to T — 1, 0.5, 0.2, and 0.05 in order of 
thinnest to thickest. 

peak of the triple resonance {Q = coi). Note further 
that the full response is not just an enhancement of the 
nonresonant features (which are apparent when the in- 
cident photon frequency becomes large), but that the 
shape of the response can change dramatically due to 
resonant effects. This is most apparent when the in- 
cident photon energy is close to U. In Fig. 2, we plot 
the resonant profile at fixed Q = 0.5 as a function of 
the incident photon frequency Wi. Note how the low 
energy features change their resonant behavior from 
being centered around uJi « 0.5 — Q (dominated by 
triple-resonance effects) at high temperature to being 
centered around Wi « t/ at lower temperatures. This 
resonance of a low-energy feature, due to a higher- 
energy photon energy has been often seen in Raman 
scattering in strongly-correlated materials. 
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